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ABSTRACT 


Experiments  were  conducted  to  determine  whether  energy-density 
relations  can  be  used  to  predict  elastic-plastic  stresses  and  strains  near  a 
circular  notch  for  7075-T6  aluminum  alloys  and  ARALL  4  composites.  The 
loading  conditions  were  tension  and  four-point  bending.  Glinka  and  Neuber 
have  developed  relations  that  predict  local  inelastic  strain  response  based  on 
the  stress-strain  solution  for  small  plastic  zone  sizes.  It  has  been  shown  that 
these  relations  are  appropriate  for  simple  tension  and  in-plane  bending, 
where  stress  and  strain  are  uniform  through  the  thickness.  This  study 
investigates  the  application  of  the  Glinka  and  Neuber  relations  to  samples 
where  stress/strain  is  not  constant  though  the  thickness.  Nonconstant 
stresses /strains  are  the  result  of  out-of -plane  bending  and  laminate 
characteristics. 
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1.  INTRODUCTION 


An  important  aspect  of  fracture  analysis  and  fatigue  life  prediction  of  a 
notched  sample  is  the  relationship  between  the  load,  the  nominal  elastic 
stress,  and  the  elastic-plastic  notch  tip  strain.  Material  failure  is  often  initiated 
by  the  localized  yielding  from  the  stress  concentration  at  a  notch  root. 
Because  solutions  for  local  stresses  and  strains  in  the  inelastic  region  are 
difficult  to  obtain,  researchers  have  attempted  to  solve  for  the  inelastic 
stresses  and  strains  based  on  simpler  elastic  solutions.  A  popular  technique 
currently  used  to  obtain  accurate  solutions  is  the  finite  element  method 
(FEM).  The  drawbacks  to  this  method  are  difficulties  in  using  nonlinear 
elements  and  lengthy  computational  times,  especially  for  a  large  number  of 
cyclic  loads.  Approximate  methods  have  been  developed  that  generate  good 
results  in  selected  cases. 

This  study  considers  the  energy-density  method  of  approximating 
inelastic  notch  tip  stresses  and  strains.  The  foundation  of  this  method  is  that 
the  energy  density  in  the  elastic  range  is  nearly  the  same  as  the  energy  density 
in  the  inelastic  range.  Previous  authors  have  applied  this  method  to  several 
types  of  steels  and  strain-hardening  materials  and  to  several  types  of  notches 
(e.g.,  circular,  elliptical,  V,  keyhole,  etc.).  They  also  considered  tension,  in¬ 
plane  bending,  and  cyclic  loading.  The  stress  and  strain  in  these  tests  were 
uniform  through  the  sample  thickness. 
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In  this  study,  experimental  testing  was  conducted  on  7075-T6 
aluminum  alloy  and  aramid  aluminum  laminate  (ARALL  4)  samples  with  a 
stress  concentration  around  a  circular  notch.  By  applying  tension  to  the 
ARALL  4  sample  and  out-of-plane  bending  to  both  samples,  a  condition  of 
nonconstant  stress  and  strain  through  the  thickness  was  created.  To  analyze 
these  conditions,  several  assumptions  were  made: 

•  All  samples  were  in  a  plane  stress  condition. 

•  For  small  scale  localized  plastic  yielding  around  a  notch,  the  inelastic 
energy  density  was  the  same  as  the  elastic  energy  density. 

•  Sample  loading  was  monotonic. 

•  A  imiaxial  stress  state  was  considered  at  the  vicinity  of  the  notch  tip 
because  of  the  small  local  stress  oriented  transverse  to  loading 
direction. 

The  primary  focus  of  this  work  is  an  experimental  study  to  determine 
whether  energy-density  based  approximation  methods  can  be  applied  to  7075- 
T6  aluminum  alloys  for  out-of-plane  bending  and  ARALL  4  composites  for 
both  tension  and  out-of-plane  bending.  The  material  properties  for  each 
material  will  be  verified  and/or  established  for  this  study.  For  each  loading 
t5^e,  the  experimental  setup  will  be  described  and  the  theory  used  to  predict 
inelastic  stresses /strains  will  precede  the  results  of  the  tests. 
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IL  LITERATURE  SURVEY 

A.  IMPORTANT  ADVANCES  ON  STRESSES  AROUND  CIRCULAR 

NOTCHES 

In  1929,  Howland  [1]  used  a  successive  approximation  method  to  solve 
for  the  stresses  in  a  semi-infinite  plate,  with  a  circular  hole  mid-way  between 
the  parallel  edges,  under  tension.  He  developed  general  formulas  in  an 
infinite  series  form  that  were  used  to  plot  the  stress  variation  around  the 
hole.  Each  term  of  the  series  was  based  on  the  ratio  of  hole  diameter  to  plate 
width  and  constant  coefficients  were  derived  by  Howland. 

In  1939,  Dumont  [2]  tested  a  large  plate  containing  an  open  circular  hole 
by  subjecting  the  plate  to  tmiform  bending  normal  to  the  plane  of  the  plate. 
He  found  that  stress  concentrations  occurred  at  the  edge  of  the  hole.  Further, 
stress  became  constant  at  a  distance  of  2.5  times  the  hole's  diameter  from  the 
edge  of  the  hole. 

In  1952,  Hill  [3]  stated  that  in  practice  the  radius  of  curvature  of  the 
notch  root  must  be  two  or  three  times  the  plate  thickness  if  the  theory  of 
generalized  plane  stress  is  to  be  a  satisfactory  approximation. 

In  1974,  Peterson  [4]  published  a  book  of  stress  concentration  factors  for 
numerous  plate  geometries.  Related  to  this  study  was  the  plot  of  stress 
concentration  factors  for  the  transverse  bending  of  a  finite-width  plate 
containing  a  circular  hole.  He  did  not  plot  the  stress  gradient  ahead  of  a  hole. 
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but  rather  the  stress  concentration  factor  at  the  edge  of  the  hole  for  ratios  of 
hole  diameter  to  plate  width  from  X  =  0.0  to  0.2. 

In  1990,  Prasad  and  Shuart  [5]  developed  a  closed  form  solution  for  the 
moment  distribution  around  holes.  Their  analytical  study  considered 
symmetric  composite  laminates  subjected  to  bending  moments. 

B.  HISTORY  OF  ENERGY  DENSITY  METHOD 

The  most  frequently  used  relation  for  predicting  local  elastic-plastic 
stresses  and  strains  near  notches  was  developed  by  Neuber  [6]  in  1961.  More 
recent  studies  have  shown  that  Neuber's  rule  can  lead  to  considerable  over¬ 
estimations  of  the  elastic-plastic  strains  at  the  notch  tip. 

In  1964,  Theocaris  and  Marketos  [7]  used  a  birefringent  coating  method 
complemented  with  the  electrical  analogy  method  to  determine  the 
distribution  of  principle  stresses  and  strains  around  a  circular  notch.  They 
considered  AL57S,  a  strain-hardening  aluminum  alloy,  for  monotonic 
loading  up  to  generalized  plastic  flow. 

In  1968,  Hutchinson  [8]  showed  that  the  strain  energy  in  the  plastic 
zone  ahead  of  a  crack  tip  is  the  same  as  that  calculated  based  on  the  linear 
elastic  stress-strain  analysis  for  material  characterized  by  a  bilinear  stress- 
strain  behavior.  In  other  words,  the  energy  density  at  the  notch  root  does  not 
change  significantly  if  the  localized  plasticity  is  surroimded  by  predominately 
elastic  material. 

In  1981,  Molski  and  Glinka  [9]  used  Hutchinson's  energy  density  theory 
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and  made  the  assumption  that  the  theory  could  be  applied  to  notches  and 
material  exhibiting  nonlinear  stress-strain  behavior.  A  uniaxial  stress  state 
was  considered  at  the  notch  tip  subject  to  cyclic  loading. 

In  1985,  Glinka  [10]  applied  the  energy  density  method  for  a  uniaxial 
stress  state  at  the  notch  tip  undergoing  monotonic,  nonproportional  loading. 
Theoretical  strain  approximations  were  compared  to  experimental  strains 
previously  measured  by  Theocaris  and  Marketos  [7].  Glinka  compared 
inelastic  notch  tip  strains  for  several  notches  and  materials.  This  study  uses 
the  uniaxial  stress  state  with  monotonic  loading  and  applies  it  to  tension 
loading  of  an  ARALL  4  composite  and  out-of-plane  bending  of  7075-T6 
aluminum  and  ARALL  4  composite. 

In  Glinka's  [11]  second  paper  of  1985,  he  modified  the  strain  energy 
density  method  by  adding  a  factor  Cp  that  accounts  for  the  increase  in  plastic 

zone  size  due  to  local  yielding.  This  factor  becomes  important  when  strains 
exceed  1.0  percent. 

Also  in  1985,  Hoffmann  and  Seeger  [12]  developed  a  method  of 
estimating  elastic-plastic  notch  stresses  and  strains  for  a  multiaxial  stress  state 
in  the  vicinity  of  the  notch  tip.  This  method  applied  to  both  monotonic  and 
cyclic  loading. 

In  1988,  Glinka,  Ott  and  Nowack  [13]  used  the  equivalent  strain  energy 
density  method  to  determine  the  elastic-plastic  stresses  and  strains  in  notches 
under  a  plane  strain  condition.  Results  were  compared  with  FEM  solutions. 
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In  1992,  Sharpe,  Yang  and  Tregoning  [14]  conducted  an  in-depth  review 
of  the  Neuber  and  Glinka  relations  for  monotonic  loading.  They  determined 
that  in  general,  the  Neuber  method  was  best  suited  for  a  plane  stress 
condition  because  it  was  conservative.  The  Glinka  plane  stress  method 
provided  an  upper  bound  of  stress.  For  large  plastic  deformation  regions,  the 
Glinka  and  Neuber  methods  bounded  the  experimental  results. 

In  1995,  Lee,  Chang  and  Wong  [15]  improved  Hoffmann's  and  Seeger's 
method  [12]  for  monotonic  loading  by  taking  a  general  approach.  Multiaxial 
notch  stresses  and  strains  were  estimated  based  on  elastic  notch  stress 
solutions  for  in-phase  or  out-of-phase  loading.  The  authors  also  reported  that 
in  general,  the  Glinka  and  Neuber  methods  underestimate  the  notch  strain  in 
the  high  yield  deformation  region  and  overestimate  the  local  strain  in  the 
low  plasticity  zone. 
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III.  MATERIAL  PROPERTY  VERIFICATION 


This  chapter  describes  the  preliminary  tests  and  results  that  were 
conducted  on  two  types  of  material:  7075-T6  aluminum  alloy  and  Aramid 
Aluminum  Laminate  (ARALL  4).  Preliminary  tests  were  required  to  verify 
and  establish  material  properties  for  the  material  provided  because  that 
information  was  required  for  tests  described  in  subsequent  chapters. 

A.  7075-T6  ALUMINUM  ALLOY 

1.  Material 

The  aluminum  alloy  used  in  this  study  was  cut  from  one  sheet  of 
7075-T6  aluminum.  This  aluminum  sheet  metal  was  also  identified  by  the 
Government  Mill  Form  and  Condition  designation  as  QQ-A-250-12.  Initial 
dimensions  of  the  large  sheet  were:  54  in.  (137.2  cm)  x  46  in.  (116.8  cm)  x  0.125 
in.  (0.318  cm). 

2.  Experimental  Apparatus 

Tension  tests  were  conducted  on  four  test  coupons  to  verify  the  7075-T6 
aluminum  sheet  metal  material  properties  used  in  this  study  as  compared 
with  published  7075-T6  aluminum  material  properties.  The  test  coupon 
design  and  the  testing  methods  used  were  governed  by  American  Society  for 
Testing  of  Materials  (ASTM)  standard  test  methods  [16].  Figure  1  shows  the 
dimensions  of  the  tension  test  coupons  used  for  material  verification.  The 
methods  listed  in  [16]  were  used  to  determine  the  material  properties  of  yield 
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strength  (offset=0.2%)  and  ultimate  tensile  strength.  The  testing  methods 
used  to  determine  Young's  modulus  and  Poisson's  ratio  were  listed  in  [17] 
and  [18],  respectively. 


A  -  length  of  gage  section 

3  1/2  in. 

B  -  length  of  grip  section 

2  in. 

C  -  width  of  grip  section 

3/4  in. 

G  -  gage  length  for  extensiometer 

3  1/2  in. 

L  -  overall  length 

8  in. 

T  -  thickness 

1/8  in. 

W  -  width 

1/2  in. 

Figure  1.  Tension  Test  Coupon  Dimensions. 


Test  coupon  width  and  thickness  were  measured  using  a  Mitutoyo 
Digital  Micrometer,  accurate  to  0.00005  in.  (0.001  mm).  The  90  degree  strain 
gage  rosettes  used  were  Measurements  Group  Inc.  type  CEA-13-125UT,  120 
ohm,  gage  factor  2.085  +  0.5%  (longitudinal)  and  2.115  +  0.5%  (transverse). 
The  Instron  dynamic  extensiometer  used  had  a  3.5  in.  (8.9  cm)  gage  length 
and  a  ±  0.2  in.  (0.51  cm)  range.  The  Instron  Model  4507  tension/ compression 
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test  machine  was  used  with  a  200  KN  load  cell.  Load  and  displacement  data 
were  obtained  using  the  Instron  Model  4500  data  acquisition  tower  with  the 
Instron  series  IX  automated  materials  testing  software,  version  5.28.  Strain 
gage  readings  were  acquired  through  the  Measurements  Group  P-3500  strain 
indicator /SB-100  switch  and  balance  unit  combination.  These  two  units  were 
used  to  acquire  all  the  strain  gage  readings  in  following  tests.  The  equipment 
set  up  used  for  this  and  subsequent  tests  is  pictured  in  Figure  2.  All  data 
recorded  from  the  tests  was  manually  entered  into  MATLAB  programs  that 
generated  the  graphs  presented  in  this  study. 

3.  Results 

Table  1  shows  the  comparison  of  the  experimentally  determined 
material  properties  with  published  material  properties  [19]  for  bare  7075-T6 
aluminum  alloy.  The  experimental  results  provided  in  Table  1  were  based  on 
data  and  calculations  from  three  of  the  four  test  coupons.  Three  test  coupons 
were  cut  in  a  longitudinal  orientation  from  the  aluminum  sheet.  Because 
transverse  mechanical  properties  of  many  products,  particularly  tensile 
strengths  and  ductility  in  the  short  transverse  direction,  can  be  less  than  those 
in  the  longitudinal  direction  [20],  one  test  coupon  was  cut  from  the  sheet  in  a 
transverse  direction  to  verify  isotropic  behavior  in  the  aluminum  sheet.  A 
comparison  of  transverse  and  longitudinal  material  properties  is  listed  in 
Table  2.  Because  transverse  and  longitudinal  material  properties  were  in 
good  comparison,  the  aluminum  sheet  was  considered  isotropic. 
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Experimental 

Published 

Percent 

Difference 

MPa 

ksi 

MPa 

ksi 

Yield 

Strength 

(offset=0.2%) 

505 

73.3 

503 

73.0 

0.4 

Ultimate 

Tensile 

Strength 

576 

83.6 

572 

83.0 

0.7 

Voting's 

Modulus 

69,600 

10,100 

71,000 

10,300 

2.0 

Poisson's 

Ratio 

0.37 

0.33 

12.1 

Table  1.  Comparison  of  Experimental  and  Published  Material  Properties  for 

7075-T6  Aluminum  Alloy. 


Sample  1 
(transverse) 

Sample  2 
(longitudinal) 

Sample  3 
(longitudinal) 

Sample  4 
(longitudinal) 

Yield  Strength 
(offeet=0.2%) 
(ksi) 

69.9 

73.2 

73.4 

73.3 

Ultimate 

Tensile 

Strength 

(ksi) 

83.2 

83.5 

83.7 

83.7 

Voting's 

Modtdus 

(ksi) 

9,928 

10,249 

9,944 

10,169 

Poisson's  Ratio 

0.373 

0.373 

0.369 

0.366 

Table  2.  Comparison  of  Transversly  Oriented  Sample  1  to  Longitudinally 

Oriented  Samples  1,  3,  and  4. 
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a.  Young's  Modulus 

For  Young's  modulus  calculations,  each  test  coupon  was  tension 
tested  three  times  in  the  elastic  range  by  cyclic  loading,  as  per  Reference  21. 
The  modulus  values  calculated  for  each  of  the  three  runs  were  averaged, 
resulting  in  an  average  modulus  value  for  that  test  coupon.  The  Young's 
modulus  used  in  this  study  was  the  average  modulus  value  of  the  three  test 
coupons  oriented  in  the  longitudinal  direction. 

During  the  tensile  tests,  strain  data  was  collected  from  two 
sources.  One  source  was  a  strain  gage  rosette  mounted  on  the  geometric 
center  of  the  test  coupon.  The  other  source  was  an  extensiometer.  Yotmg's 
modulii  were  calculated  using  stress  and  both  strain  data  sources.  As  before, 
modulus  values  calculated  for  each  rim  were  averaged  to  result  in  a  modulus 
value  for  each  test  coupon.  Averaging  the  longitudinally  oriented  test 
coupons'  modulii  based  on  extensiometer  data  resulted  in  Young's  modulus 
of  10,269  ksi  (70.6  GPa).  The  average  longitudinal  Young's  modulus  based  on 
strain  gage  data  was  10,121  ksi  (69.6  GPa),  a  1.5%  difference.  Reference  22 
describes  the  standard  method  of  recording  strain  data  with  the  use  of 
extensiometers.  In  order  to  prevent  damage  to  the  extensiometers  by 
subjecting  them  to  impact  forces  that  occur  at  ultimate  failure,  they  were 
removed  just  prior  to  sample  failure.  The  only  means  of  obtaining  strain 
data  up  to  ultimate  failure  was  with  the  use  of  strain  gages. 
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b.  Poisson's  Ratio 

Reference  23  provided  the  standard  method  of  experimentally 
determining  Poisson's  ratio  using  extensiometers.  Since  the  use  of  strain 
gages  and  extensiometers  produced  similar  results  for  Young's  modulus, 
strain  gage  based  data  was  also  used  for  determining  Poisson's  ratio.  The 
Instron  data  acquisition  system  and  automated  testing  software  computed  the 
ultimate  tensile  strength.  Yield  strengths  (0.2%  offset)  were  taken  from  stress- 
percent  elongation  graphs. 

c.  Nonlinear  Coefficients 

A  relation  most  often  used  for  strain  in  the  nonlinear  region  is 
the  Ramberg-Osgood  relation  [24] 

a  a  -i  ,  , 

e=£+(Y)"  (3.1) 

where  K  is  the  strength  coefficient  and  n  is  the  strain-hardening  index.  K  and 
n  were  experimentally  determined  from  the  results  of  the  tension  tests  used 
for  material  property  verification.  For  most  metals,  a  log-log  plot  of  true 
stress  versus  true  plastic  strain  is  modeled  as  a  straight  line  [25].  However,  the 
relationship  often  deviates  from  linearity  at  low  strains  and/or  high  strains 
[26].  The  slope  of  the  line  is  n  and  the  intercept  with  the  log  true  stress  axis  is 
K.  Figure  3  shows  the  plot  of  log  true  stress  versus  log  true  plastic  strain  for 
sample  2.  Similar  plots  were  seen  with  samples  3  and  4.  K  and  n  values  were 
calculated  for  samples  2,  3,  and  4  and  averaged.  It  is  obvious  that  the  data  is 
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not  linear  throughout  the  plastic  range,  but  it  is  satisfactory  to  calculate  K  and 
n  for  the  strain  range  over  which  the  log-log  plot  is  linear.  Therein  lies  the 
question  of  which  region  of  the  plot  should  be  considered  linear.  Two  cases 
of  linear  approximations  were  investigated.  Case  1  considers  the  data  points 
of  log  true  plastic  strain  between  -4  and  -2.4.  Case  2  considers  the  data  points 
of  log  true  plastic  strain  between  -3.2  and  -2.4.  A  least  squares  approximation 
of  the  data  for  both  cases  is  shown  in  Figure  3.  The  results  were  K  =  110.13  ksi, 
n  =  0.0790  for  case  1  and  K  =  114.88  ksi,  n  =  0.0946  for  case  2.  It  will  be  shown 
later  that  the  case  2  linear  approximation  for  K  and  n  will  provide  slightly 
better  results.  Thus,  the  average  K  and  n  values  used  in  this  study  are  114.88 
ksi  and  0.0946,  respectively. 

4.  Statistical  Analysis 

A  statistical  study  was  employed  to  determine  the  errors  in  the  least 
squares  method  of  approximation  to  Young's  modulus  calculations. 
Correlation  coefficients  (r)  were  calculated  for  each  run  and  averaged  for  each 
test  coupon.  The  test  coupons'  averaged  coefficients  were  then  averaged  to 
obtain  a  single  correlation  coefficient.  Correlation  coefficients  indicate  the 
reliability  of  the  curve  fit  to  the  data  and  a  perfect  curve  fit  has  a  coefficient 
equal  to  1.0.  The  average  correlation  coefficient  for  longitudinally  oriented 
stress-strain  plots  in  the  elastic  region  was  0.99979.  When  the  data  is  basically 
linear,  as  was  the  case  in  these  tensile  tests,  the  correlation  coefficient  is  not  a 
very  sensitive  indicator  of  the  precision  of  the  data.  McClintock  found  that 
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Figure  3.  Log  of  True  Plastic  Stress-Stram  Behavior  for  Sample  2.  Strength 


Coefficient  K  and  Strain  Hardening  Index  n  were  determined  using  Case  2. 


.5 


J-  2 

(l-r^)^  is  a  better  indicator  of  the  curve  fit  quality  [27].  Here,  (1-r  = 

0.0146.  This  value  indicates  that  the  vertical  standard  deviation  of  the  data  is 
only  1.46  percent  of  the  total  vertical  variation  caused  by  the  straight-line 
relationship  between  stress  and  strain  [28]. 

B.  ARAMID  ALUMINUM  LAMINATE 

1.  Material 

An  aramid  aluminum  laminate  (ARALL)  has  also  been  considered  in 
this  study.  ARALL  is  classified  as  a  polymer  matrix  composite  material.  The 
specific  fiber-reinforced  composite  used  was  ARALL  4,  which  has  a  5/4 
configuration.  There  were  five  layers  of  2024-T8  aluminum  alloy  and  four 
layers  of  unidirectional  and  continuous  Kevlar  49  aramid  fibers.  The 
laminate  was  symmetric  and  balanced  with  respect  to  the  test  direction.  The 
2024-T8  aluminum  lamina  were  isotropic  and  all  the  aramid  laminas  were 
oriented  in  the  same  direction.  Two  sheets  of  the  composite,  measuring  12 
inches  by  12  inches  each,  were  provided  for  testing. 

2.  Experimental  Apparatus 

Four  samples  were  initially  cut  from  the  sheets  provided.  Two 
samples  were  machine  finished  to  same  dimensions  as  the  7075-T6 
aluminum  samples  used  for  the  material  verification  tensile  tests  (see  Figure 
1).  One  of  these  samples  had  the  fibers  oriented  parallel  to  the  loading 
direction  and  the  other  had  the  fibers  oriented  transverse  to  the  loading 
direction.  This  design  allowed  for  the  calculation  of  longitudinal  and 
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transverse  modulii.  The  other  two  samples  were  machine  finished  to  a 
rectangular  shape  of  1  in.  (2.54  cm)  by  9  in.  (22.86  cm).  These  samples  had 
fibers  oriented  only  in  the  longitudinal  direction  and  were  used  in  flexural 
tests  to  determine  the  bending  modulus. 

The  equipment  used  for  the  laminate  material  property  tensile  tests 
was  the  same  as  that  used  for  the  7075-T6  aluminum  material  property  tests. 
Figure  4  shows  an  ARALL  sample  undergoing  a  tension  test.  For  the  flexural 
tests,  a  four-point  bending  apparatus  was  constructed  to  determine  bending 
modulus  and  it  is  shown  in  Figure  5.  Four-point  loading  at  the  quarter  points 
was  the  load  configuration,  which  is  a  typical  loading  condition  for  high- 
modulus  materials.  The  support  span  was  5  5/8  in.  (143  mm)  and  the  load 
span  was  2  13/16  in.  (71.5  mm).  The  sample  length  was  7  1/4  in.  (188  mm) 
and  laminate  thickness  was  3/32  in.  (2.4  mm).  Sample  length  was  based  on 
an  ASTM  recommended  60:1  span-to-thickness  ratio  for  high-modulus 
composites.  Deflection  was  measured  using  a  Laboratory  Devices  Co.  dial 
gage,  accurate  to  0.001  in.  The  strain  gages  used  were  Measurements  Group 
Inc.  type  CEA-13-250UN,  350  ohm,  gage  factor  2.12  +  0.5%.  For  both  tensile 
and  flexural  tests,  the  Instron  test  machine  was  used  in  the  manual  load 
control  mode  to  record  load,  strain  and  deflection  data. 
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Figure  4.  Tension  Test  of  an  ARALL  4  Sample. 
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Figure  5.  Four-Point  Bending  Apparatus  for  Flexural  Tests. 

3.  Results 

The  reference  used  for  the  tensile  test  to  determine  Young's  modulus 
and  Poisson's  ratio  was  ASTM  3039  [29]  and  the  reference  used  for  the  flexural 
test  to  determine  the  bending  modulus  was  ASTM  790  [30]. 
a.  Young's  Modulus 

Two  samples  were  tested  to  determine  a  smeared  Young's 
modulus.  This  was  the  modulus  for  the  laminate  rather  than  a  modulus  for 
each  of  the  different  lamina.  One  sample  was  used  to  calculate  (Young's 

modulus  in  the  fiber  direction)  and  the  other  was  used  to  calculate  E2 
(Young's  modulus  transverse  to  the  fiber  direction).  Two  strain  gages  were 
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placed  on  each  sample.  One  gage  was  oriented  in  the  fiber  direction  and  the 
other  gage  was  oriented  transverse  to  the  fiber  direction.  Although  each 
sample  was  tested  to  ultimate  failure,  only  data  in  the  linear  elastic  range  was 
used  to  calculate  modulus  values.  For  the  sample  with  fibers  oriented 
longitudinal  to  the  loading  direction.  Figure  6  shows  the  stress-strain  curves 
in  the  longitudinal  and  transverse  directions.  Figure  7  shows  the  stress-strain 
curves  in  the  longitudinal  and  transverse  directions  for  the  case  of  fibers 
oriented  transverse  to  the  loading  direction.  Figure  8  shows  the  plot  of  stress 
versus  strain  for  longitudinal  and  transverse  directions  in  the  linear  elastic 
range.  A  MATLAB  program  was  written  to  solve  for  the  moduli!  in  a  least 
squares  sense.  The  best  fit  approximation  of  the  slope  of  the  data  in  the  elastic 
region  for  Figures  6  and  7  was  the  modulus  value.  For  the  longitudinal 
modulus,  Ej  =  9,100  ksi  (62.6  GPa)  and  for  the  transverse  modulus,  =  6,600 

ksi  (45.4  GPa).  These  modulii  show  that  a  sample  with  fibers  oriented 
longitudinal  to  loading  direction  is  stiffer  than  a  sample  with  fibers  oriented 
transverse  to  loading  direction. 

A  comparison  of  the  calculated  E^  value  to  published  elastic 

modulii  (tension)  for  2024  aluminum  and  Kevlar  49  fibers  is  made.  Because 
elastic  modulii  are  dependent  on  laminate  thickness  and  lamina 
orientation,a  standard  elastic  modulus  for  the  ARALL  4  material  provided 
was  not  available.  The  published  elastic  modulus  (tension)  for  2024 
aluminum  is  10,500  ksi  (72.4  GPa)  [31]  and  for  Kevlar  49  fibers  is  18,100  ksi  (125 
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Stress  [ksi] 
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Figure  7.  Stress-Strain  Behavior  in  Longitudinal  and  Transverse  Directions 
for  Fibers  Oriented  Transverse  to  Loading  Direcion. 
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Figure  8.  Longitudinal  and  Transverse  Modulii  in  the  Elastic  Range. 


GPa)  [32].  Although  the  elastic  modulus  for  Kevlar  49  is  high,  the  lamina  did 
not  consist  of  a  100  percent  volume  fraction  of  Kevlar  49  fibers.  These  lamina 
contained  an  undetermined  amount  of  resin,  interface  bonding  material,  and 
volume  fraction  of  Kevlar  49  fibers  much  less  than  100  percent.  The  low 
volume  fraction  of  Kevlar  49  fibers  effectively  lowers  the  laminar  elastic 
modulus  to  a  percentage  equal  to  the  volume  fraction.  Thus,  laminate  elastic 
modulus  values  greatly  depend  on  the  volume  fraction  of  the  fibers  within 
the  lamina.  It  follows  that  =  9,100  ksi  (62.6  GPa)  is  a  reasonable  laminate 

elastic  modulus. 

b.  Bending  Modulus 

An  alternative  approach  to  tensile  testing  of  unidirectional 
pol5mieric-matrix  composites  involves  the  bending  of  beam-type  samples  [33]. 
The  flexural  test  method  was  used  to  determine  a  smeared  bending  modulus 
for  the  overall  laminate  vice  individual  2024  aluminum  or  aramid 
aluminum  lamina  modulii. 

Two  samples  were  also  used  in  flexural  tests  to  determine  the 
bending  modulus.  Eg.  The  purpose  of  using  two  samples  was  to  achieve 

repeatability  of  the  data.  Since  the  results  for  each  sample  were  nearly 
identical,  only  one  sample's  data  will  be  presented.  Figure  9  shows  the  plot  of 
load  versus  deflection  in  the  elastic  range.  For  the  test.  Eg  =  11,200  ksi  (77.0 

GPa).  The  bending  modulus  was  calculated  using 
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(3.2) 
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where  L  =  length  of  support  span  [in.  (mm)],  b  =  with  of  beam  [in.  (mm)],  t  = 
thickness  of  beam  [in.  (mm)],  and  m  =  slope  of  load-deflection  curve  in  elastic 
range  [Ibf/in.  (N/mm)]. 

For  laminate  composites,  the  tensile  modulus  and  bending 
modulus  are  expected  to  be  slightly  different.  With  highly  anisotropic 
laminates,  the  bending  modulus  is  strongly  dependent  on  the  lamina¬ 
stacking  sequence  [34].  A  tensile  modulus  is  not  stacking-sequence 
dependent.  In  the  flexural  test  conducted,  the  bending  modulus  Eg  was  19 

percent  higher  than  the  longitudinal  modulus  E^.  This  can  be  explained  by 

the  different  stiffnesses  of  the  lamina.  The  2024  aluminum  lamina  are  stiffer 
than  the  aramid  fiber  lamina.  The  outermost  lamina  of  the  composite  were 
the  2024  aluminum  lamina.  The  stiffer  outermost  lamina  results  in  higher 
bending  stiffness  but  has  no  effect  on  the  tensile  stiffness, 
c.  Poisson's  Ratio 

The  strain  data  taken  from  each  sample  used  in  the  Young's 
modulus  section  was  used  to  determine  Poisson's  ratio.  Two  Poisson's  ratios 

were  calculated,  corresponds  to  the  sample  with  fibers  oriented 

longitudinally  to  the  loading  direction,  while  is  for  the  sample  with  fibers 
oriented  transverse  to  the  loading  direction.  Figure  10  shows  Poisson's  ratio. 
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^12'  versus  load  in  the  elastic  region.  The  average  value  is  drawn 

bythe  horizontal  line  and  calculated  by 


v 


avg 


n 


(3.3) 


where  Ejj  is  the  longitudinal  strain,  e^.  is  the  transverse  strain,  and  n  is  the 
total  number  of  data  points  in  the  elastic  region.  This  calculation  was  used  to 
find  =  0.323.  Reference  31  reports  Poisson's  ratio  as  0.33  for  2024 

aluminum.  This  compares  nicely  with  the  calculations  for  Since  the 

laminae  containing  longitudinally  oriented  Kevlar  49  fibers  have  almost  no 
transverse  stiffness,  the  laminate  Poisson's  ratio  was  essentially  due  to  the 
effects  of  the  2024  aluminum  lamina. 

The  relation 


E, 


(3.4) 


can  be  used  to  predict  since  the  other  values  are  known.  The  calculation 
yields  v^^  =  0.234.  For  comparison,  the  plot  of  experimentally  determined 
Poisson's  ratio,  v^^,  versus  load  in  the  elastic  region  is  shown  in  Figure  11. 
Averaging  the  ratios  results  in  v^^  =  0.270,  which  is  15  percent  higher  than  the 
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predicted  value.  This  difference  can  be  attributed  several  factors.  One  factor  is 
that  the  experimental  precision  used  to  determine  Young's  modulus  is 
related  to  gage  factor  errors  and  load  readout  errors.  Another  factor  is  that  the 
plot  of  v^j  values  in  Figure  11  shows  a  wide  range  of  data  point  scatter. 

However,  the  predicted  v^^  value  does  fall  within  the  range  of  data  point 

scatter.  The  factors  combine  to  give  a  slightly  higher  v^^  than  expected.  As 

expected.  Figures  10  and  11  also  show  that  Poisson's  ratio  is  higher  for  fibers 
oriented  in  the  loading  direction  than  in  the  transverse  direction. 
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Poisson  Ratio 
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Figure  11.  Poisson's  Ratio  in  the  Elastic  Range  for  Fibers  Oriented  Transverse 

to  Loading  Direction. 
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IV.  TENSION  TEST 


Tension  tests  on  circular  notched  samples  were  conducted  for  both  the 
7075-T6  aluminum  and  the  aramid  aluminum  laminate  samples.  The 
purpose  of  the  tests  on  the  7075-T6  aluminum  was  to  verify  Glinka's  previous 
work  in  Reference  10  and  to  validate  the  equipment  used  in  this  study. 
Tension  testing  on  the  ARALL  4  was  done  to  determine  if  the  energy-density 
method  of  calculating  inelastic  stress-strain  near  a  circular  notch,  imder  plane 
stress  conditions,  could  be  applied  to  an  ARALL  4  composite.  This  chapter 
discusses  the  equipment  apparatus  used  for  tension  testing,  the  theory  of  the 
energy-density  method  based  on  two  different  equations,  the  method  of 
computations  for  each  material,  and  the  results  of  the  comparison  of  Glinka's 
and  Neuber's  theoretical  inelastic  stress-strain  approximations  to 
experimental  stresses  and  strains.  Test  results  indicated  that  the  Glinka  and 
Neuber  methods  worked  as  well  as  expected  for  the  7075-T6  aluminum 
sample  but  not  for  the  ARALL  4  sample.  Good  theoretical  inelastic  strain 
predictions  can  be  made  using  modified  Glinka  and  modified  Neuber 
methods. 

A.  EXPERIMENTAL  APPARATUS 

Figure  12  shows  the  geometry  and  dimensions  of  the  circular  notched 
samples  used  for  the  tension  tests.  The  length  of  the  7075-T6  aluminum 
sample  was  15  in.  and  the  length  of  the  ARALL  4  sample  was  12  in.  Both 
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Figure  12.  Circular  Notched  Sample  Dimensions  and  Strain  Gage  Locations. 
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sample  types  had  machine  finished  edges.  The  circular  notch  was  milled 
through  in  the  7075-T6  aluminum  sample  and  drilled  through  in  the  ARALL 
sample.  The  circular  notch  was  located  in  the  geometric  center  of  the  sample. 
Also  shown  in  Figure  12  is  the  location  and  orientation  of  the  strain  gages. 
The  10-element  strain  gages  used  were  Measurements  Group  Inc.  type  CEA- 
09-020PF,  120  ohm,  gage  factor  2.00  +  3.0%.  The  single  element  strain  gages 
used  were  Measurements  Group  Inc.  type  CEA-13-032UW,  120  ohm,  gage 
factor  2.11  +  1.0%.  The  two  single  element  strain  gages  were  positioned  back- 
to-back  only  on  the  7075-T6  aluminum  sample.  Having  the  two  single  gages 
positioned  back-to-back  and  midway  between  the  outer  edge  of  the  hole  and 
edge  of  the  sample  permitted  the  comparison  of  strains  on  both  sides  of  the 
sample.  Greater  strains  on  one  side  of  the  sample  would  indicate  the 
presence  of  bending  stresses  caused  by  misaligned  sample  gripping.  Tension 
tests  were  conducted  on  the  Material  Test  System  (MTS)  Model  810 
tension/ compression  test  machine.  The  load  cell  cartridge  was  rated  at  25 
kips  and  the  maximum  deflection  cartridge  was  rated  at  1  inch.  The  MTS  test 
machine  was  selected  for  testing  because  the  grip  width  was  2  in.,  which  was 
wide  enough  to  extend  across  the  entire  width  of  the  circular  notched 
samples.  A  circular  notched  sample  undergoing  a  tension  test  is  pictured  in 
Figure  13.  For  the  7075-T6  aluminum  sample,  the  length  of  the  sample 
gripped  by  the  MTS  machine  was  2  in.  on  each  end.  The  ARALL  4  sample 
gripped  length  was  1.25  in.  on  each  end.  The  MTS  machine  was  manually 
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Figure  13.  Tension  Test  of  a  Circular  Notched  Sample. 
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load  controlled.  Using  a  10,000  sec/in.  loading  speed,  the  loading  was 
periodically  paused  to  allow  recording  of  load  and  strain  data.  Samples  were 
loaded  to  ultimate  failure. 

B.  THEORY  OF  ENERGY-DENSITY  METHOD 

1.  Based  on  the  Ramberg-Osgood  Relation 

In  Reference  10,  Glinka  described  the  energy-density  method  for 
calculating  inelastic  stresses-strains  near  notches.  To  assist  the  reader,  the 
energy-density  method  of  calculating  inelastic  stress-strain  near  a  circular 
notch  under  plane  stress  conditions  is  presented  below. 

To  calculate  the  energy  density,  a  stress-strain  relation  must  be  applied. 
Well  known  and  often  used,  the  Ramberg-Osgood  relation  is 

e=  e,+ Ep  =  -|  +  (^)~.  (4.1) 

The  strain  energy  density  distribution  ahead  of  a  notch  tip,  W^,  can  be 
calculated  as 

e. 

W,  =  js„de„  (4.2) 

0 

where  S..  is  the  local  elastic  stress  tensor  and  e..  is  the  local  elastic  strain 
1] 

tensor.  For  localized  small-scale  plastic  yielding,  the  energy  density  in  the 
plastic  zone  is  nearly  equal  to  that  in  the  elastic  zone.  Thus,  the  energy 
density  in  the  plastic  zone  is  equal  to  that  calculated  on  the  basis  of  an 

elastic  solution 
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(4.3) 


e  £. 

1  ) 


W,  =  or  JSijdejj  =  JCijdEij 

0  0 


where  a.,  is  the  local  inelastic-plastic  stress  tensor  and  £..  is  the  local  elastic- 


plastic  strain  tensor. 

The  elastic  stress  at  the  notch  tip  is 

S,  =  K,S„ 


(4.4) 


where  is  the  local  elastic  stress  component  in  the  y  direction,  is  the  stress 
concentration  factor,  and  is  the  nominal  stress.  For  uniaxial  stress  at  the 
notch  tip,  equation  (4.3)  can  be  written  in  the  form 

(4.5) 


(K,  ■ 


2E 


y 

2E  ■ 


This  states  that  elastic  strain  energy  density  at  the  notch  tip  is  equal 
to  the  product  of  the  strain  energy  density  due  to  the  nominal  stress 
and  the  square  of  the  stress  concentration  factor 

K?W3„  =  W„.  (4.6) 


With  localized  yielding  at  the  notch  tip,  the  energy  density  should  be 
calculated  with  respect  to  the  nonlinear  stress-strain  relation  (4.1),  then 
equation  (4.6)  takes  the  form 


K 


n  +  1 


(4.7) 


When  the  nominal  stress  has  exceeded  the  proportional  limit,  equation 
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(4.6)  takes  the  form 


but  the  right  hand  side  does  not  represent  energy  density,  is  the  Neuber 

determined  inelastic  stress.  The  difference  between  Neuber's  rule  and 
Glinka's  energy-density  relation  (4.8)  is  the  before  the  second  term  in  each 

side  of  (4.8)  is  replaced  by  ay  in  (4.10). 
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2.  Based  on  the  Bilinear  Relation 

It  will  be  shown  later  in  this  chapter  that  the  predicted  theoretical 
inelastic  notch  strains  for  the  ARALL  4  composite  are  unsatisfactory  when  the 
Glinka  and  Neuber  methods  are  used.  The  Glinka  and  Neuber  equations  can 
be  modified  to  account  for  the  bilinear  stress-strain  behavior  of  the  ARALL  4 
composite. 

In  the  upper  linear  region  of  the  bilinear  stress-strain  diagram,  (4.7) 
becomes 


K 


2 

t 


2E  2E, 


for  (a  <a) 
y  1 


(4.11) 


where  is  the  modulus  in  the  upper  linear  region  [ksi]  and  Oj  is  the  stress  at 

the  intersection  of  each  linear  region  [ksi].  Theoretical  inelastic  notch  strains 
for  the  modified  Glinka  method  are  approximated  using 


ay-Oi 


(4.12) 


where  Cj  is  the  strain  at  the  intersection  of  each  linear  region. 

Neuber's  method  may  be  modified  to  accovmt  for  the  bilinear  stress- 
strain  behavior,  such  that 


r^2 


(4.13) 


where  represents  the  stress  from  Neuber's  method  [ksi].  Theoretical 
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inelastic  notch  strains  for  the  modified  Neuber's  method  are  also  foimd  using 
(4.12). 

C  METHOD  OF  COMPUTATIONS 

To  compare  experimental  stress-strain  with  theoretical  stress-strain 
determined  by  the  energy-density  method,  the  experimental  data  was  entered 
into  a  MATLAB  program.  The  different  material  types  required  different 
equations  to  solve  for  inelastic  notch  strains. 

1.  7075-T6  Aluminum  Alloy 

For  nominal  stresses  lower  than  the  yield  stress,  theoretical  elastic 
strains  were  calculated  using  Hooke's  Law.  Equations  (4.8)  and  (4.10)  were 
used  to  compute  Glinka's  and  Neuber's  theoretical  inelastic  stresses  based  on 
the  nominal  stresses  calculated  in  (4.9).  Because  (4.8)  and  (4.10)  are  nonlinear, 
the  bisection  numerical  technique  was  employed  to  solve  for  the  theoretical 
inelastic  stresses.  Theoretical  inelastic  strains  were  calculated  by  substituting 
the  theoretical  inelastic  stresses  from  (4.8)  and  (4.10)  into  (4.1). 

2.  ARALL  4  Composite 

In  the  elastic  region,  strains  were  calculated  using  Hooke's  Law  and 
Glinka's  equation  (4.7).  Inelastic  strains  were  calculated  using  Glinka's 
equation  (4.10).  These  equations  produced  unsatisfactory  results.  Better 
results  were  obtained  when  Glinka's  and  Neuber's  methods  were  modified  in 
(4.11)  and  (4.13).  Inelastic  stresses  were  calculated  from  (4.11)  and  (4.13)  and 
substituted  into  (4.12)  to  obtain  inelastic  strains. 
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D.  RESULTS 


1.  7075-T6  Aluminum  Alloy 

a.  Theoretical  Circular  Notch  Strain 
Figure  14  shows  that  the  Ramberg-Osgood  relation  produces  a  good 
stress-strain  approximation  to  experimental  data  for  7075-T6  aluminum. 
However,  the  decision  as  to  which  data  points  to  consider  to  determine  the 
strength  coefficient  K  and  strain  hardening  index  n  can  affect  the  results. 
Figure  15  shows  the  comparison  of  Glinka's  theoretical  notch  strain  to 
experimental  notch  strain  for  tensile  loading.  The  experimentally 
determined  stress-strain  curve  is  added  to  Figure  15  for  reference.  Two 
theoretical  approximations,  cases  1  and  2,  are  shown.  In  Case  1,  K  =  110.13  ksi 
and  n  =  0.0790.  In  Case  2,  K  =  114.88  ksi  and  n  =  0.0946.  As  previously 
mentioned  in  Chapter  III  and  now  apparent,  case  2  K  and  n  material  property 
values  provide  a  closer  approximation  to  experimental  results.  During  elastic 
behavior,  the  experimental,  theoretical,  and  stress-strain  curves  are  the  same. 
The  deviation  of  the  experimental  and  theoretical  curves  from  the  stress- 
strain  curve  occurs  at  the  beginning  of  inelastic  behavior.  In  the  inelastic 
range,  the  results  confirm  Glinka's  results  [10]  in  that,  "The  difference 
between  calculated  and  measured  notch  strains  were  smaller  than  10  percent 
...".  For  the  7075-T6  aluminum  sample,  the  maximum  calculated  and 
measured  strain  difference  was  11.8  percent  in  the  inelastic  region.  This 
experiment  confirmed  the  Lee  et  al.  [15]  claim  that  Glinka  s  rule 
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Figure  14.  Ramberg-Osgood  Stress-Strain  Approximation  to  Experimental 
Data  for  7075-T6  Aluminum  Alloy. 
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Figure  15.  Comparison  of  Theoretical  and  Experimental  Circular  Notch 
Strains  in  7075-T6  Aluminum  Alloy  for  Two  Different  Cases  of  K  and  n 

Values. 


underestimates  the  notch  strain  in  the  high  yield  deformation  region  and 
overestimates  the  local  strain  in  the  low  plasticity  zone. 

Figure  16  shows  the  comparison  of  Glinka  and  Neuber  methods 
for  determining  theoretical  circular  notch  strain  in  7075-T6  aluminum.  As 
Glinka  [10]  indicated,  Neuber's  rule  underestimates  notch  strain.  However, 
Sharpe  et  al.  [14]  reports  that  the  Neuber  rule  is  the  best  single  model  for 
predicting  notch  strains  in  a  plane  stress  condition.  This  is  because  the 
Neuber  rule  is  more  conservative  than  the  Glinka  rule.  For  large  plastic 
deformation  regions,  the  Glinka  plane  stress  model  gives  an  upper  bound  of 
strain.  The  sample  reached  ultimate  failure  at  a  tensile  load  of  13,270  lbs;  the 
corresponding  nominal  stress  was  85.0  ksi. 

b.  Stress  Concentration  Factors 

The  plot  of  theoretical  and  experimental  circular  notch  strains 
proved  to  be  very  sensitive  to  the  selection  of  the  elastic  stress  concentration 
factor  Kj.  The  elastic  stress  concentration  factor  gradient  ahead  of  the  notch 

was  calculated  using  Howland's  method  of  coefficients  [1].  The  ratio  of  hole 
diameter  to  plate  width  was  X  =  0.334.  Howland  only  calculated  coefficients 
for  the  ratios  of  =  0.1,  0.2,  0.3,  0.4  and  0.5.  Thus,  the  coefficients  used  in  this 

study  had  to  be  interpolated  from  Howland's  table  of  coefficients.  The  stress 
concentration  factor  gradient  away  from  the  notch  is  shown  in  Figure  17.  The 
gradient  is  steepest  near  the  notch  and  decreases  away  from  the  notch. 
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Figure  16.  Comparison  of  Glinka  and  Neuber  Methods  for  Determining 
Theoretical  Circular  Notch  Strain  in  7075-T6  Aluminum  Alloy. 


Superimposed  on  the  gradient  curve  are  the  locations  of  individual  strain 
gages  (indicated  by  horizontal  error  bars)  from  the  10-element  strip  gage.  The 
average  stress  concentration  factor  over  the  strain  gage  range  is  indicated  by  a 
circle.  Next  to  the  hole,  =  2.26  using  Howland's  method.  Integrating  the 

gradient  over  the  range  of  the  first  strain  gage's  measuring  range  to  find  the 
average,  =  2.04.  The  stress  concentration  factor  used  for  the  7075-T6 

aluminum  tension  test  results  was  K^=  1.90.  Although  this  fell  within  the 

range  of  the  first  strain  gage,  it  was  slightly  lower  than  the  average  because  the 
first  gage  was  located  a  small  distance  away  from  the  edge  of  the  notch. 

Figure  18  shows  the  theoretical  strain  gradient  ahead  of  the 
circular  notch  for  two  nominal  stresses  in  the  elastic  range:  =  34.0  ksi  and 

=  19.2  ksi.  In  [10],  Glinka's  application  of  the  energy-density  method 

resulted  in  an  underestimation  of  theoretical  inelastic  strains  in  strain 
gradient  plots.  This  study  found  that  theoretical  inelastic  strains  based  on 
Howland's  solution  were  slightly  higher  than  experimental  strains.  Figure  19 
shows  strain  gradients  ahead  of  the  circular  notch  for  several  nominal 
stresses.  At  low  nominal  stresses,  the  gradient  is  nearly  linear.  However,  as 
stress  is  increased,  the  gradient  next  to  the  edge  of  the  notch  greatly  increases. 
2.  ARALL  4  Composite 

Both  Glinka's  equation  (4.8)  and  Neuber's  equation  (4.10),  that  are  used 
to  predict  theoretical  inelastic  stress,  depend  on  the  strength  coefficient,  K, 
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Figure  18.  Theoretical  and  Experimental  Elastic  Range  Strain  Gradients 
Ahead  of  the  Circular  Notch  in  7075-T6  Aluminum  Alloy. 
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and  the  strain  hardening  index,  n.  Assuming  the  stress-strain  behavior  for 
the  ARALL  4  composite  can  be  represented  by  the  Ramberg-Osgood  equation 
(4.1),  K  and  n  must  be  calculated  and  applied  to  (4.8)  and  (4.10).  It  is  evident 
from  Figure  6  that  the  stress-strain  behavior  of  the  ARALL  4  composite  is 
bilinear.  Figure  20  shows  the  plot  of  log  true  plastic  stress  versus  log  true 
plastic  strain.  The  method  for  calculating  K  and  n  for  the  ARALL  4  composite 
was  similar  to  that  for  the  7075-T6  aluminum.  The  straight  line  in  the  upper 
linear  region  of  the  curve  was  used  to  calculate  K  =  889  ksi  and  n  =  0.570. 

Figure  21  shows  the  comparison  of  experimental  and  theoretical  notch 
strains  using  Glinka's  method.  Glinka's  method  greatly  overestimated  the 
inelastic  notch  strains.  Neuber's  method  was  not  shown  because  it  produced 
nearly  the  same  overestimation  for  inelastic  strains  that  Glinka's  method  did. 
Labeled  A-B-C  is  a  strain  discontinuity  for  the  transition  from  the  linear 
elastic  region  to  the  inelastic  region.  Inelastic  strains  were  calculated  using 
(4.8).  If  elastic  strains  are  calculated  using  Hooke's  Law,  the  jump  in  strain 
going  from  elastic  to  inelastic  behavior  is  shown  as  the  line  A-C.  Using  (4.7) 
in  the  elastic  region,  the  corresponding  jump  in  strain  is  shown  as  line  B-C. 
The  strain  discontinuity  between  (4.7)  and  (4.8)  can  be  attributed  to  the 
equations'  dependence  on  K  and  n. 

Figure  22  shows  the  Ramberg-Osgood  stress-strain  approximation  to 
experimental  data  for  the  ARALL  4  composite.  Clearly,  the  Ramberg-Osgood 
equation  does  not  model  bilinear  stress-strain  behavior.  Consequently,  the 


49 


Glinka  equations  produced  unsatisfactory  results.  The  stress  concentration 
factor  used  in  the  ARALL  tension  test  results  was  2.10.  The  sample  reached 
ultimate  failure  at  a  tensile  load  of  8,550  lbs;  the  corresponding  nominal  stress 
was  72.8  ksi. 

The  strain  gradients  ahead  of  the  notch  for  several  nominal  stresses  are 
shown  in  Figure  23.  The  behavior  of  the  gradients  as  the  stress  was  increased 
was,  as  expected,  similar  to  that  of  the  7075-T6  aluminum. 

The  prediction  for  theoretical  inelastic  notch  strains  can  be  greatly 
improved  by  using  the  modified  Glinka  (4.11)  and  modified  Neuber  (4.13) 
equations.  The  comparison  of  the  modified  Glinka  and  modified  Neuber 
methods  for  the  ARALL  4  composite  in  tension  is  shown  in  Figure  24.  These 
two  new  methods  give  a  bound  to  the  experimental  data.  The 
overestimation  in  inelastic  strains  using  the  modified  Neuber  method 
provides  a  more  conservative  approximation.  A  slight  underestimation  in 
notch  strains  resulted  from  the  modified  Glinka  method. 
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Figure  21.  Comparison  of  Experimental  and  Theoretical  Notch  Strains  Using 
Glinka's  Method  for  an  ARALL  4  Composite  in  Tension. 
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Figure  22.  Ramberg-Osgood  Stress-Strain  Approximation  to  Experimental 

Data  for  ARALL  4  Composite. 
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Figure  24.  Comparision  of  Modified  Glinka  and  Modified  Neuber  Methods 

for  ARALL  4  Composite  in  Tension. 
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V.  FLEXURAL  TEST 


Flexural  tests  were  conducted  on  both  7075-T6  aluminum  and  ARALL 
4  samples.  Four-point  loading  was  applied  to  the  quarter  points  of  the 
samples.  This  created  out-of-plane  bending  stresses  around  the  circular 
notch.  Flexural  testing  was  conducted  to  determine  whether  the  Glinka  and 
Neuber  relations  established  for  tensile  loading  could  be  used  predict  local 
inelastic  strains  for  bending.  This  chapter  discusses  the  experimental 
apparatus  used  in  conducting  the  tests,  the  flexural  test  procedure,  the  theory 
used  to  determine  inelastic  bending  stresses  for  both  materials,  and  the  test 
results  for  the  7075-T6  aluminum  and  ARALL  4  composite. 

A.  EXPERIMENTAL  APPARATUS 

The  design  of  the  circular  notched  samples  used  in  the  flexural  tests 
were  the  same  as  that  used  in  the  tension  tests.  Sample  sizes  and  dimensions 
listed  in  Figure  12  also  apply  to  flexural  test  samples.  The  same  10-element 
strip  gages  were  used  but  the  single  back-to-back  strain  gages  were  omitted. 
Also,  the  Instron  test  machine  used  for  material  property  verification  was 
used  for  flexural  tests.  However,  in  order  to  take  incremental  load  and  strain 
data,  the  Instron  was  operated  in  the  manual  load  control  mode.  A  four- 
point  bend  test  apparatus  was  manufactured  in  accordance  with  ASTM 
specifications  [30]  and  the  dimensions  are  shown  in  Figure  22.  The 
equipment  setup  is  pictured  in  Figure  23. 
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A  -  Load  Span  Length 

7075-T6  A1  -  7  15/16  in. 

ARALL4-  5  1/16  in. 

B  -  Support  Span  Length 

7075-T6  A1  -  12  1/16  in. 

ARALL  4  -  10  1/8  in. 


Figure  25.  Dimensions  of  Four-Point  Bend  Test  Apparatus. 
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Figure  26.  Flexural  Test  of  Circular  Notched  Sample  Subjected  to 

Four-Point  Loading. 
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B.  BEND  TEST  PROCEDURE 


A  sample  was  placed  on  the  base  of  the  beam  test  apparatus  and 
symmetrically  aligned.  Each  strain  gage  of  the  10-element  strip  gage  was 
zeroed.  The  top  of  the  beam  test  apparatus,  which  weighed  9.6  lbs,  was  placed 
on  the  sample  and  symmetrically  aligned.  This  condition  accounted  for  the 
first  data  point.  Incremental  loads  were  then  applied  to  the  sample  to  obtain 
the  remaining  data  points.  For  the  7075-T6  aluminum  sample,  the  test  was 
concluded  when  the  sample  deflected  approx.  2.1  inches,  almost  touching  the 
base  of  the  beam  test  apparatus.  At  this  point,  the  corresponding  load  was  365 
lbs.  The  maximum  load  reached  before  the  ARALL  4  sample  deflected  into 
the  base  of  the  bend  test  apparatus  was  410  lbs.  Because  the  ARALL  4  sample 
were  shorter  in  length,  it  could  withstand  a  greater  load. 

C  THEORY 

As  tensile  loading  is  increased  for  a  circular  notched  sample,  plastic 
yielding  first  occurs  at  the  edge  of  the  notch  and  progresses  away  from  the 
notch,  towards  the  edge  of  the  sample.  At  any  location  away  from  the  notch, 
the  amount  of  plastic  yielding  at  that  distance  will  be  the  same  through  the 
sample  thickness.  Just  as  in  tensile  loading,  the  plastic  yielding  for  bending 
begins  at  the  edge  of  the  circular  notch  and  progresses  outward.  However,  at  a 
given  location  away  from  the  notch,  as  loading  is  increased  plastic  yielding 
begins  on  the  surface  of  the  sample  and  progresses  inward,  through  the 
thickness,  towards  the  neutral  axis.  A  small  region  of  plastic  yielding  is  often 
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referred  to  as  localized  plastic  yielding.  For  small  localized  plastic  regions, 
where  nominal  bending  stress  exceeds  yield  stress  a^,  the  energy  density 


at  the  notch  tip  can  be  calculated  based  on  bending  elastic  stress. 

1.  7075-T6  Aluminum  Alloy 

The  same  Glinka  relation  presented  in  Chapter  IV  for  localized 
yielding  at  the  notch  tip  under  tensile  loading  can  be  used  for  bending 


-L 

n+1  K  ^  ■ 


The  Neuber  relation  for  bending  is  similar  to  (4.7)  but  the  in  the 
second  term  on  the  right  hand  side  is  replaced  with  y,  giving 


‘  2E  2E  2  K 


(5.1) 


Nominal  bending  stress  is  calculated  by 

6M 

S„  =  ^ - 

t^p-2d) 


(5.2) 


where  M  is  the  bending  moment  [Ibf-in],  t  is  the  sample  thickness  [in.],  D  is 
the  plate  width  [in.]  and  d  is  the  notch  diameter  [in.]. 

2.  ARALL  4  Composite 

The  modified  Glinka  and  modified  Neuber  equations  presented  in 
Chapter  IV  can  be  also  used  for  bending.  The  modified  Glinka  equation  is 
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The  modified  Neuber  equation  is 
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D.  RESULTS 

1.  7075-T6  Aluminum  Alloy 

Figure  27  shows  the  comparison  of  Glinka  and  Neuber  methods  for 
determining  theoretical  inelastic  notch  strain  in  out-of-plane  bending.  The 
Glinka  approximation  was  much  closer  to  the  experimental  strains  than  the 
Neuber  approximation.  For  low  and  moderate  local  inelastic  strains,  the 
Glinka  method  was  an  overestimation  but  in  the  high  yield  region,  the 
Glinka  method  tended  to  underestimate  the  local  inelastic  strain. 

Figure  28  shows  the  load  versus  deflection  curve  for  7075-T6 
aluminum  in  out-of-plane  bending.  The  tangent  to  the  linear  region  of  the 
curve  was  shown  to  identify  the  departure  from  linearity.  Figure  29  shows 
several  strain  gradients  ahead  of  the  circular  notch.  These  gradients  are 
slightly  different  when  compared  to  those  for  the  7075-T6  aluminum  samples 
in  tension  loading.  The  slope  of  the  gradient  next  to  the  notch  was  lower  for 
bending  than  tension.  As  in  the  tension  test,  the  strain  gradients  in  bending 
approached  a  constant  value  away  from  the  notch. 

The  stress  concentration  factor  used  in  Figures  27  and  29  was  1.20. 
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Figure  27.  Comparison  of  Glinka  and  Neuber  Methods  for  Determining 
Theoretical  and  Experimental  Circular  Notch  Strain  for  7075-T6  Aluminum 

in  Out-of-Phase  Bending. 


Strain 


Figure  29.  Strain  Gradient  Ahead  of  the  Circular  Notch  for  Several  Nominal 
Stresses  for  7075-T6  Aluminum  in  Out-of-Plane  Bending. 
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Peterson  has  plotted  several  curves  for  stress  concentration  factors  subject  to 
transverse  bending  of  a  finite- width  plate  with  a  circular  hole  [35].  The 
experimental  circular  notch  diameter  (a)  to  sample  width  (W)  ratio  was  0.34. 
Peterson  plotted  the  stress  concentration  factor  curves  for  from  0  to  0.2.  If 

the  appropriate  curve  is  extrapolated  to  -^  =  0.34,  the  resulting  stress 

concentration  factor  next  to  the  notch  is  approximately  1.40.  As  discussed  in 
Chapter  IV,  the  strain  gage  was  located  small  distance  from  the  edge  of  the 
notch.  The  strain  gage  averaged  the  stress  concentration  factor  gradient  over 
the  gage's  measured  area.  The  15  percent  difference  in  experimental  and 
theoretical  stress  concentration  factors  can  also  be  attributed  to  a  slight  off- 
centering  of  the  circular  notch  by  0.03  in.  (0.76  mm)  during  the  machining 
process. 

2.  ARALL  4  Composite 

Figure  30  shows  the  comparison  of  Glinka  and  Neuber  methods  for 
determining  theoretical  and  experimental  circular  notch  strain  for  ARALL  4 
composite  in  out-of-plane  bending.  Glinka's  equation  (4.7)  was  used  within 
the  elastic-plastic  region  to  predict  theoretical  notch  strains.  This  was  justified 
by  Glinka's  statement  in  Reference  10,  "...  within  some  range  of  the  nominal 

bending  stress  >  Cy,  only  localized  plastic  yielding  can  occur  in  the  analyzed 

section,  whereas  under  tensile  loading  plastic  yielding  of  the  whole  section 
would  occur."  The  discontinuity  line  A-B  shows  the  transition  from  elastic 
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strains  calculated  using  Hooke's  Law  and  inelastic  strains  calculated  using 
Neuber's  method.  Although  the  Glinka  method  achieved  a  closer 
approximation,  both  methods  greatly  overestimated  local  inelastic  strains.  At 
low  inelastic  yielding,  the  two  methods  predicted  a  strain  response  that 
exceeded  the  material  stress-strain  diagram.  As  in  Chapter  IV,  the  equations' 
dependence  on  the  strength  coefficient  (K)  and  strain  hardening  index  (n) 
from  the  Ramberg-Osgood  relation  resulted  in  unsatisfactory  results. 

Figure  31  shows  the  comparison  of  the  modified  Glinka  and  the 
modified  Neuber  methods  for  out-of-plane  bending.  These  methods  greatly 
improved  the  inelastic  strain  approximation  as  compared  to  Figure  30.  Both 
modified  methods  underestimated  experimental  notch  strain  at  large  plastic 
deformations.  For  small  plastic  deformations,  the  modified  Neuber  method 
was  a  slightly  conservative  approximation. 

In  Figure  32,  the  load-deflection  curve  for  the  ARALL  4  composite  in 
out-of-plane  bending  is  shown.  Figure  33  shows  several  strain  gradients 
ahead  of  the  circular  notch  for  the  ARALL  4  composite  in  out-of-plane 
bending.  The  stress  concentration  factor  used  in  the  ARALL  4  tests  was  1.05. 
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Figure  30.  Comparison  of  Gliirka  and  Neuber  Methods  for  Determining 
Theoretical  and  Experimental  Notch  Strain  for  ARALL  4  Composite  in  Out- 

of-Plane  Bending. 
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Figure  31.  Comparison  of  Modified  Glinka  and  Modified  Neuber  Methods  for 
ARALL  4  Composite  in  Out-of-Plane  Bending. 
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Bending. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


Because  the  behavior  of  inelastic  stresses /strains  around  a  circular 
notch  is  complex,  an  accurate  solution  has  not  been  found.  Approximate 
methods  have  been  proposed  that  provide  reasonable  results.  Specific 
findings  of  this  study  include: 

•  Neuber's  method  always  predicted  inelastic  circular  notch  strains 
that  were  larger  than  those  using  Glinka's  method. 

•  As  expected,  Glinka's  method  of  approximating  theoretical  inelastic 
notch  strain  produced  good  results  for  the  7075-T6  aluminum  in  tension. 
Theoretical  inelastic  strains  were  within  15  percent  of  the  experimental  data. 

•  The  Glinka  method  provided  good  theoretical  strain  predictions  for 
a  7075-T6  aluminum  alloy  subjected  to  out-of-plane  bending. 

•  Both  the  Glinka  and  Neuber  methods  demonstrated  unsatisfactory 
results  when  applied  to  an  ARALL  4  composite  subjected  to  tension  and  out- 
of-plane  bending.  Because  separate  equations  were  used  to  predict  strains  in 
the  elastic  and  inelastic  region,  a  strain  discontinuity  was  created  in  the 
transition  from  linear  elastic  to  inelastic  behavior.  ARALL  4  has  a  bilinear 
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stress-strain  curve  and  Glinka's  and  Neuber's  method  are  based  on  the 


Ramberg-Osgood  relation.  This  relation  does  not  represent  the  composite's 
bilinear  stress-strain  behavior  well. 

•  The  modified  Glinka  and  modified  Neuber  methods  produced  good 
results  for  predicting  theoretical  inelastic  notch  strains  in  an  ARALL  4 
composite  subjected  to  out-of-plane  bending. 

•  For  an  ARALL  4  composite  in  tension,  the  modified  Glinka  and 
modified  Neuber  methods  can  be  used  to  bound  experimental  notch  strains. 
Using  the  standard  Glinka  and  Neuber  methods,  bounding  of  experimental 
notch  strains  was  also  observed  for  the  7075-T6  alumimun  sample  in  tension. 

•  Theoretical  notch  strain  predictions  for  samples  in  tension  were 
better  than  samples  in  out-of-plane  bending.  This  was  true  for  the  standard 
and  modified  Glinka  and  Neuber  methods. 

•  Previous  studies  found  Neuber's  method  produced  conservative 
predictions  for  theoretical  inelastic  notch  strains.  For  the  ARALL  4  composite 
in  out-of-plane  bending,  the  modified  Neuber  method  predicted  inelastic 
notch  strains  that  were  non-conservative. 
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Future  research  in  the  inelastic  stress  distribution  around  notches 
should  consider  the  following  recommendations: 

•  Due  to  the  large  stress  concentration  factor  gradient  next  to  a  circular 
notch,  a  precise  stress  concentration  factor  at  the  notch  tip  is  difficult  to 
obtain.  Precision  could  be  improved  by  using  larger  sample  widths  and  notch 
diameters  with  small  strain  gages.  However,  the  problems  associated  with 
increased  sample  lengths  and  grip  configuration  at  the  ends  for  tension  tests 
requires  careful  planning.  To  record  large  inelastic  strains,  the  strain  gage  grid 
must  be  as  large  as  possible  to  avoid  internal  gage  failure.  The  trade  off  is  that 
strain  gages  tend  to  average  the  strain  over  the  area  covered  by  the  grid.  For  a 
nonuniform  strain  distribution,  the  reported  strain  measurement  will  less 
than  the  maximum  strain.  Consequently,  the  precise  determination  of  the 
stress  concentration  factor  at  the  notch  tip  will  be  more  difficult. 

•  Instead  of  using  strain  gages.  Moire  interferometry  could  be  used  to 
measure  the  local  displacement/strain  more  accurately. 

•  A  finite  element  analysis  could  be  conducted.  FEA  would  provide 
highly  accurate  results  and  incorporate  a  degree  of  flexibility.  Different 
sample  and  notch  geometries  could  be  modeled  for  different  sample  loadings. 
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